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Non existence of a phase transition for the Penetrable Square Wells in one dimension2 

1. Introduction 

The Penetrable Square Well (PSW) model in one dimension was first introduced in 
[1] as a good candidate to describe star polymers in regimes of good and moderate 
solvent under dilute conditions. The issue of Ruelle' s thermodynamic stability was 
analyzed and the region of the phase diagram for a well defined thermodynamic limit 
of the model was identified. A detailed analysis of its structural and thermodynamical 
properties where then carried through at low temperatures [2] and high temperatures. 

m 

The problem of assessing the existence of phase transitions for this one dimensional 
model had never been answered in a definitive way. Several attempt to find a gas-liquid 
phase transition were carried through using the Gibbs Ensemble Monte Carlo (GEMC) 
technique [H El El [71 E] but all gave negative results. Now it is well known that in 
three dimensions the Square Well (SW) model admits for a particular choice of the well 
parameters a gas-liquid transition. [9] As the van Hove's theorem shows, [lOl [121 1131 tH] 
this disappears in one dimension. Nonetheless the PSW model in one dimension, being 
a non nearest neighbors fiuid, is not analytically solvable and since we have no hard 
core the van Hove's theorem does not hold anymore. It is then interesting to answer the 
question whether a phase transition is possible for it. We should also mention that we 
also used the GEMC technique to probe for the transition in the three dimensional PSW 
and we generally found that for a given well width there is a penetrability threshold 
above which the gas-liquid transition disappears. 

In the present work we use the Kastner, Schreiber, and Schnetz (KSS) theorem 
p^ [15] to give strong analytic evidence for the absence of any phase transition for this 
fiuid model. 

The argument hinges on a theorem of Szego [16] on Toeplitz matrices and can be 
applied to a large class of one dimensional fiuid models and complement the van Hove's 
theorem. 

The paper is organized as follows: in Section [2] we state the KSS theorem for the 
exclusion of phase transitions, in Section [3] we describe the PSW model, in Section [H 
we show numerically that the PSW model satisfies KSS theorem, in Section [5] we show 
analytically that the PSW model satisfies the KSS theorem, the conclusive remarks are 
presented in Section [61 

2. The KSS theorem 

The Kastner, Schreiber, and Schnetz (KSS) theorem [HI [15] states the following. 



Theorem KSS: Let Vn '■ ^n ^ ^'^ — t- M be a smooth potential; an analytic 
mapping from the configuration space Fjv onto the reals. Let us indicate with 'H^(q) 
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the Hessian of the potential Indicating with qc the critical points (or saddle points) of 
VAr(q) (i.e. VqVAr|q=qc = Oj, with fc(qc) their index (the number of negative eigenvalues 
of 'H'^ (cic) ) . Assume that the potential is a Morse function (i.e. the determinant of the 
Hessian calculated on all its critical points is non zero). Whenever F^v is noncompact, 
assume Vn to be "confining", i.e. liniA^oo VAr(Aq) = 00, VO 7^ q G T^. Consider the 
Jacobian densities, 



where 



and 
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{^c\[Vn{c[c)/N = v]A [fc(q,) = /(mod4)]} . (3) 

Then a phase transition in the thermodynamic limit is excluded at any potential en- 
ergy in the interval {v ~ e,v + e) if: (i.) the total number of critical points is limited 
by exp{CN), with C a positive constant, (ii.) for all sufficiently small e the Jacobian 
densities are jiiy) < +00 for I = 0, 1, 2, 3. 



Generally the number of critical points of the potential grows exponentially with 
the number of degrees of freedom of the system. The fact that the total number of 
critical points is limited by an exponential is thought to be generically valid. [17] We 
then assume that for Morse potentials the first hypothesis of the theorem is satisfied. 
So the key hypothesis of the theorem is the second one, which can be reformulated as 
follows: for all sequences of critical points qc such that limAr_>.oo VAr(qc)/A^ = v, we have 

qJI^^O. (4) 



lim I det Ti 

Af-s>oo 



A^, 



3. The PSW model 



The pair potential of the PSW model can be found as the / — >■ 00 limit of the following 
continuous potential 

0;(r) = a[b - tanh(/(r - 1))] + c[tanh(/(r - A)) + 1] , (5) 

where a = (e^ + ea)/2, b = (e^ — ea)/(er + ea), c = ea/2, with e^ a positive constant which 
represent the degree of penetrability of the particles, €„ a positive constant representing 
the depth of the attractive well, and A = 1 + A, with A the width of the attractive 
square well. The Penetrable Spheres (PS) in one dimension are obtained as the A — )■ 
limit of the PSW model. In the hmit of e^ — )■ 00 the PSW reduces to the SW model. 
The PSW model is Ruelle stable for e^/e^ > 2(n + 1) with n<A<n + l. [Il[3] 



Non existence of a phase transition for the Penetrable Square Wells in one dimension^ 



Let us consider a pair potential of the following form 

2 




1 — cos I 27r— 
Ij 



(6) 



Note that this pair potential is periodic of period L and flat at the origin, $;(0) = 0. 
Moreover in the large L limit $i(r) k. (j)i(r'^). In Fig. [1] we show this potential for 
different choices of the smoothing parameter /. 
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Figure 1. Shows the potential $;(|a;|) for L ':^ 1. In the plot we used e^ = 5,ea 
1, A = 4, and L = 10, at two values of the smoothing parameter I. 



4. Absence of a phase transition 

In this section we will apply the KSS theorem to give numerical evidence that there is 
no phase transition for the PSW model introduced above. 
The total potential energy is 

1 ^ 
Vn{cl) = -Y,M\^.-x,\), (7) 

where q = (xi, a;2, . . . , xn)- If liniAr_j.oo VAr(q)/A^ = v one finds er/2 — ea < v < +oo. 

The saddle points qs = (x^, x|, . . . , x%) for the total potential energy (VqVAf = 0), 
can be various. We will only consider critical point of the following kind: equally spaced 
points at fixed density p = N/L, 

Vp, ^ = 0,l,2,...,iV-l. (8) 



x- 



Here we can reach 



lim lV(qp)/iV 

W— s>oo 



-■p 5 



(9) 
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where for large A^ and up to an additive constant — 0/(O)/2 we have, 

cos(^)]) . (10) 

If p ^ 1, in the big A^ hniit we can approximate the sum by an integral so that 

^ ' dx . 11 



vr Jo Vl-(l-x)2 
keeping in mind that L = N/p and A^ is big we find in the / — )■ oo limit 

Vp ~ — {er[- arcsin(l - z)]]^^^^^^^""^'^ - £,[- arcsin(l - z)]^^^^^^^',]} 
« 2p[e,- 6,(71 -l)] = t;°, (12) 

where we used for small z, arcsin(l — z) = n/2 — \plz + 0[z'^/^]. 
For small p in the I — )■ oo limit you get, 

Vp = e,/2 , p < l/V\ (13) 

Vp = e^/2 - ta , l/VX < p < 1 (14) 

For intermediate values of the density you will get a stepwise function of the density. 
A graph of Vp is shown in Fig. [2l 

Other stationary points would be the ones obtained by dividing the interval L into 
p = N/a (a > 1) equal pieces and placing a particles at each of the points x^ '^ = iL/p, 
i = 0, . . . ,p — 1. By doing so we can reach limTv-s-oo VN{(iN,p)/N = VN^p where up to an 
additive constant — 0/(O)/2 we have 







, 2TTi 
1 — COS 

p 



(15) 



We then immediately see that for p ^ a, limAr^oo "^A^.p = Vp but for small p, vn,p > Vp. 
The Hessian T-L^jiq) = d^VN^c^/dxidxj calculated on the saddle points of the first 
kind can be written as 



nUci,)= -$;'(r,,), 1^3 , (16) 



N 



where $; (r) is the second derivative of $i(r) and rjj = \i — j\/p- 

So the Hessian calculated on the saddle point is a circulant symmetric matrix with 
one zero eigenvalue due to the fact that we have translational symmetry xf = xf ± n/p 
for any i and any integer n. In order to break the symmetry we need to fix one point 
for example the one at x^. So the Hessian becomes a (A^ — 1) x (A^ — 1) symmetric 
Toeplitz matrix (non circulant anymore) which we call 'H*^^~^''(qp). 
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Figure 2. Shows the behavior of Vp as a function of the density p ioi N = 100, 200, 
and 300 when e^ = 5, Ca = 1, and A = 2 with I — 100. Also the theoretical prediction 
fp at big densities (Eq. (|12p ) is shown. Notice that at fixed TV, Vp will saturate to 
sa Ntr for 4(L/27r)2 < 1 or p > N/-k. 



In Fig. |3]we have calculated the | det 'H^(qp)|^/^ as a function of A^ at p = N/L 
fixed for ea = 1, e,. = 5, A = 1, and / = 10. One can see that the normalized determinant 
of the Hessian does not go to zero in the large N limit. So the Kastner, Schreiber, and 
Schnetz (KSS) criteria [HI [15] is not satisfied and a phase transition is excluded. The 
same holds for the PS model. 

In Fig. Hlwe show the dependence of | det7{^(qp)|^/^ on density for different choices 
of A^. 
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Figure 3. Shows the behavior of | detH^(qp)|^/^ as a function of N at two different 
densities. Here we chose Eq = 1, e^ = 5, A = 1, and I = 10. 
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Figure 4. Shows the behavior of | defH^(qp)p^ as a function of p for various N . 
Here we chose €„ = l,er = 5,A = 1, and I = 10. Notice that for p < 1/vA then 
^^(Qp) ~ ^nd also the normalized determinant is very small. While the approach 
to zero at large densities is an artifact of the finite sizes of the systems considered. 



A system where there is a phase transition has been proved to be the self- 
gravitating ring (SGR) [18] where (pscRii") = — l/^/r + 2{L/2'K)'^e. |j In this case 
one finds w^ = -p2^/2/^A{2/e), with A{x) = J^^"^ dO {1 + xsin^ 6)-^/^ They use 
Hadamard upper bound to the absolute value of a determinant to prove that indeed 
lim7v-s>oo I det'H^(qp)|^/^ = 0. In Fig. [5] we show this numerically for a particular 
choice of the parameters. Actually this result could be expected from what will be 
proven in the next section, as in the large N limit for any finite e, (psGR = o{l/N) and 



det?^^(q, 



|1/7V 



o{l/N). This is a confirmation that theorem KSS is not violated. 



5. Limit of the normalized determinant 



In this section we will give analytical evidence that there cannot be a phase transition 
for the PSW model. 

We need to apply to our case, Szego's theorem p^ for sequences of Toeplitz ma- 
trices which deals with the behavior of the eigenvalues as the order of the matrix goes 
to infinity. In particular we will be using the following Proposition. 



Proposition: Let Tn = {tkj\k-,j = 0, 1, 2, . . . ,n — 1} be a sequence of Toeplitz 



matrices with fL = t'?_- such that T = lim„^oo Tn and tk 



kj 



k-j 



lim„^ooifc fork = 0,1,2.... 



I With this choice the pair potential ^sgr. would be 27rp times the pair potential in the paper of 

Nardini and Casetti. (TB] 

§ Note that there is an error in the paper of Nardini and Casetti. [THj 



Non existence of a phase transition for the Penetrable Square Wells in one dimensions 




1e-018 



0.001 



0.010 



1/N 



Figure 5. Shows the behavior of | defH^(qp)p^ as a function of N for fixed p — 1 
in a bilogarithmic plot. Here we chose e = 0.1. 



Let us introduce 



f{x)= Y. t>^^''^ ^ xG[0,27r] 



k=—oo 






such that 



lim I detT^I^/" = lim | detTj^/'^ = exp ( — / In |/(x)| dx ] , (20) 

71— ^oo n—^OD \ 2,71 Jq 

as long as the integral of In | / (x) | exists finite. 



1 



27r 



(18) 



Then there exists a sequence of Toeplitz matrices Tn = {tkj\k,j = 0,1,2, ... ,n — 1} 
with ti^j = tk_j and 

1 r'i-K 



(19) 



If the Toeplitz matrix is Herniitian then t_k = t*^ and / is real valued. If moreover 
The Toeplitz matrix is symmetric then t_fc = tk and additionally f{x) = /(27r — x). 

By choosing T/v = 'H^(qp) and calling tf_j = 'H^(qp) we have in the iV — )■ oo limit, 
with L = N/p (p constant), tk = lim.N^oot^ and 



{N-l)/2 

f(x)= lim I 2 V t^cos(kx)+t^ 



N 



AT odd \ fc=l 



2 y^ tk cos{kx) + to J 



(21) 



fc=i 
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tf = -$;'(A;/p), A; = l,2,...,(iV-l)/2, (22) 

(N-l)/2 

t^ = -2 5^ tf , (23) 

k=l 

So that /(O) = 0. Notice that in this case the sequence of matrices 'H'^{qp) does not 
coincide with the sequence used in the Proposition, only the hmiting matrix for large 
A^ coincides. But since Szego's theorem states the limit of the normalized determinant 
exists it should be independent from the sequence chosen. An additional support to the 



Proposition is presented in [Appendix A 



Now in order to prove the absence of a phase transition we need to prove that 
Jp'^ln |/(a;)| dx does not diverge to minus infinity. That is we must control the way / 
passes through zero. In particular we do not want to have that if xo is a zero of / then 

|/(x)|~e-i/l^-^'»l", x~xo, (24) 

with a > 1, which is faster than any finite power of {x — Xq). 

Now for PSW we can write $/(r) = $™''^(r) + $*''*' (r). Choose $J"*'(r) = 
aexp(— 2Zr^) with a = (e^ + er)e^' — e^e^'^'. It is then always possible to redefine the 
starting potential ^i{r) in such a way that ^'^"'^^(r) exactly vanishes for r > Vcut > v^ 
keeping all the derivatives at r = r^ut continuous, [jjj Now in Eq. 0211) for f^°'^'^ only a Gnite 
number of k contributes to the series, namely the ones for 1 < k < pvcut- So f'^"^'^ will be 
well behaved on its zeroes. For the tail we get f*"'^\x) = —a^ixj^lx^ exp(— x^/8/). So 
that we will never have |/(a;)| going through a zero (note that the zeroes of / increase in 
number as p increases) with the asymptotically fast behavior of Eq. ( 1241) . This proves 
the absence of any phase transition for the PSW (or PS) models. 

Note that the argument continues to hold for example for the Gaussian Core 
Model (GCM) [I9] defined by 4>GCM{r) = eexp[— (r/a)^]. In this case by choosing 
0(r) = exp(— r) we get in the large L limit $(r) = exp(— r^) and the Fourier transform 
of $ (r) is —^/7Tx'^exp{—x'^/4:) which poses no problems for the zero of f{x) at x = 
(note that in this case f{x) is always positive for a; > 0). 

The argument breaks down for example if f{x) = — exp(— l/|a;|). In this case 
the pair potential will be given by $(r) ~ — /_ exp{ixr)f{x)/x'^dx, and one finds 
$(r) ~ 2[y/ —ir Ki(2\^ —ir) + \/irKi(2y/vr)], where Kn is the modified Bessel function 
of the second kind. See Fig. Elfor a plot. Also the relevant feature, in the pair potential, 
which gives the break down of the argument for the absence of a phase transition, is the 
large r behavior. Notice that in this case we numerically found out that the normalized 
determinant tend to a finite value for large A^. In accord with the fact that when the 
hypotheses of the proposition are not satisfied Eq. fl20l) looses its meaning. Considering 
the normalized determinant for the rescaled potential ^{r)/h{N), with h{N) — )■ +oo 

II Note that since the potential energy must be a Morse function (in the hypotheses of KSS theorem), 
we cannot take the tail potential $*'^*'(r) such that it exactly vanishes for r > r^ut- On the other hand 
the Gaussian decay of ^i{r) for large r is sufficient to guarantee the power law behavior of / on its 
zeroes. 
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as A^ — i- oo, we saw that it indeed tends to zero, indicating the presence of a phase 
transition. 

We simulated this model fluid and indeed we found that it undergoes a gas-liquid 
phase transition. The coexisting binodal curve is shown in Fig. [7] and in Table [1] we 
collect various properties of the two phases. We used GEMC in which two systems 
can exchange both volume and particles (the total volume V and the total number of 
particles N are fixed) in such a way to have the same pressures and chemical potentials. 
We constructed the binodal for A^ = 50 particles. In the simulation we had 2N particle 
random displacements (with a magnitude of O.Sctj, where ctj is the dimension of the 
simulation box of system i), N/10 volume changes (with a random change of magnitude 
0.1 in ln[V^i/(\^ — Vi)], where Vi is the volume of one of the two systems), and A^ particle 
swap moves. We observed that in order to obtain the binodals at different system sizes 
we had to assume a scaling of the following kind: f^N" = /SsoSO"" =constant, indicating 
that the model is not Ruelle stable (as it may be expected since it has a bounded core 
and a large attractive region), and pN = psoSO =constant, where (3^0 and pso are the 
coexistence data shown in Fig. [7] and Table [H For 50 < A^ < 100 we found a ~ 1/2, 
for A^ ^ 200 then a ^ 2/3, and for A^ ^ 300 then a ^ 3/4. 



o 




Figure 6. Shows the pair potential <l>(r) = 2 [y/~irKi {2y^—ir) + \/irKi{2\/ir)\ of the 
counterexample given in the text. We have $(0) ~ 2 and $(r) c>c sin-\/2r exp(— V^r) 
at large r. 



We then added an hard core to the potential 
<l>(r) = 



2[^/^i^Ki{2^- 



-ir] 



r < 1 

irKi(2yir)] r > 1 



(25) 



with e a positive large number, and we saw, through GEMC, that the corresponding 
fluid still admitted a gas-liquid phase transition (without A^ scaling of the densities 
p < 1) in accord with the expectation that are the large r tails of the potential that 
make this model singular from the point of view of our argument. 
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Figure 7. Shows the gas-Uquid coexistence line in the temperature density plane, 
obtained with the GEMC for A^ = 50 particles [20] interacting with the pair potential 
of Fig. El 
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0.40 


0.20 ±0.01 


1.61 ±0.03 


-0.224 ±0.009 


-0.907 ±0.007 
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-0.97 ±0.01 


0.42 


0.25 ±0.02 


1.51 ±0.02 


-0.26 ±0.01 


-0.873 ±0.008 


-0.95 ±0.01 


-0.943 ±0.008 


0.44 


0.292 ±0.007 


1.46 ±0.02 


-0.290 ±0.007 


-0.854 ±0.004 


-0.938 ±0.004 


-0.921 ±0.006 
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0.350 ±0.007 
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-0.815 ±0.006 


-0.90 ±0.01 


-0.89 ±0.02 
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0.411 ±0.007 


1.21 ±0.02 


-0.370 ±0.003 


-0.77 ±0.01 


-0.886 ±0.003 


-0.86 ±0.01 


0.50 


0.49 ±0.01 


1.04 ±0.02 


-0.420 ±0.006 


-0.71 ±0.01 


-0.87 ±0.01 


-0.862 ±0.006 



Table 1. Gas-liquid coexistence data (T, pi,Ui, fii are respectively the temperature 
the density, the internal energy per particle, and the chemical potential of the vapor 
i = V or liquid i = I phase. (3 = l/ksT and A is the de Broglie thermal wavelength.) 
from GEMC of iV = 50 particles [20]. 



For fluids with a pair potential $ given by a hard core and a — l/r° tail we can 
take the $ (r) = for r < 1 and $ (r) = —a{a — l)/r°~^ for r > 1, and the resulting / 
function (the Fourier transform of — $ ) is such that In |/(x)| has non-integrable zeros. 
So this class of models does not fall under the hypotheses pf the proposition. And it is 
well known that when 1 < a < 2 the corresponding fluid admits a phase transition |12j . 



6. Conclusions 

Using KSS theorem and a limit theorem of Szego on Toeplitz matrices we were able 
to give strong evidence for the exclusion of phase transitions in the phase diagram of 
the PSW (or PS) fluid. The argument makes use of the fact that the smoothed pair 
potential amongst the particles has an r cutoff. Even if we just considered two classes of 
stationary points, i.e. the equally spaced points and equally spaced clusters, we believe 
that our argument give strong indications of the absence of a phase transition. 

Our argument applies equally well to model fluids with large r tails in the pair 
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potential decaying in such a way that the condition of Eq. f l24p does not hold. For 
example it applies to the Gaussian Core Model. We believe this to be a rather large 
class of fluid models. 

We give an example of a model fluid which violates the condition of Eq. (!24|) and 
find through GEMC simulations that it indeed has a gas-liquid phase transition. 

Our argument does not require the fluid to be a nearest neighbor one, for which 
it is well known that the equation of state can be calculated analytically [211 |22l |23] . 
We think that our argument can be a good candidate to complement the well known 
van Hove theorem for such systems that violates the hypotheses of the hard core 
impenetrability of the particles and of the compactness of the support of the tails. 



Appendix A. Alternative support to the Szego result 

Our original matrix 'H^(qp) is a circulant matrix 
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■ h^ 1 



(A.l) 



We have numerically checked that the determinant of T-L^ {(\p) with one row and one 
column removed converges in the large N limit to the product of the non- zero eigenvalues 
of the matrix 'H^(qp). ISl 

Let us assume that A^ = 2n + l is odd. Then our matrix has the following additional 
structure 



/if =/if , ^ = l,...,n 



h 



N 



k 



N 



t = l. 



n 



"n+i n—{i—l) ' 

The eigenvalues of if ^ will be given by p4] 

N-l 



k=0 



with the additional constraint (see Eqs. ( TT6|) -( !T71) ) that 



N-l 



^0 = E ^^ = • 



fc=0 



The eigenvalues can be rewritten as follows 



^m = h^ + J2~h^e-^-' + J2h 



N 
n-(fe-l) 



g-^im(n+fc) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



fc=i 



it=i 



% We have checked numerically that this property continues to hold as long as the circulant matrix is 
a symmetric one. 
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Introducing the summation index j = n — k + 1 in the last sum we then obtain 

n 1 



n 1 

fc=l j='n 



= E tke-"^-' , (A.6) 

k=—n 

with n = {N - l)/2 and t^ = t^^ = h^ ioi k = 1,2, . . . ,n. 

We take the logarithm of the absolute value of the product of the non-zero 
eigenvalues to find 

N N 



^4^^ 



n^^- 



m=l 



N 



Now in the large A^ limit we have tk = lim7v->oo t^ for A; = 0, 1, 2, . . . and 
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i^m- Y. ^fce-t-^~/(^m) , (Ai 



V 



k=—oo 
1 ^ 

N ^-^ 

m=l 



^•t^"" 



1 

2^ 



2n 

\n\f{x)\dx, (A.9) 



'0 

where in the last passage we have transformed the sum into an integral. 
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